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ABSTRACT 


Cirrus clouds determine the radiative balance of the upper troposphere and the transport of water vapor across the tropopause. 
The representation of vertical wind velocity, W, in atmospheric models constitutes the largest source of uncertainty in the 
calculation of the cirrus formation rate. Using global atmospheric simulations with a spatial resolution of 7 km we obtain for 
the first time a direct estimate of the distribution of W at the scale relevant for cirrus formation, validated against long-term 
observations at two different ground sites. The standard deviation in W, oy, varies widely over the globe with the highest values 
resulting from orographic uplift and convection, and the lowest occurring in the Arctic. Globally about 90% of the simulated o,, 
values are below 0.1 ms~! and about one in 10* cloud formation events occur in environments with o,, > 0.8 ms~!. Combining 
our estimate with reanalysis products and an advanced cloud formation scheme results in lower homogeneous ice nucleation 
frequency than previously reported, and a decreasing average ice crystal concentration with decreasing temperature. These 
features are in agreement with observations and suggest that the correct parameterization of o,, is critical to simulate realistic 
cirrus properties. 


Introduction 


Cirrus clouds, made of ice crystals and present at low temperatures (below 235 K) and high altitudes, cover approximately 
17% of the Earth!. They control the hydrological balance of the upper troposphere and the longwave cloud feedback”. 
Outside of areas of strong convection, cirrus form by ice nucleation on aerosol particles. Human activities (e.g., transportation 
and energy production) and natural phenomena (e.g., volcanic eruptions and dust storms) that modify the concentration of 
atmospheric aerosol influence the properties of cirrus*, impacting the hydrological balance of the atmosphere and ultimately 
Earth’s climate~>-’, For example, observational studies suggest that volcanic ash leads to the increase of ice-containing clouds 
in the upper troposphere®’. There is also evidence of aerosol particles affecting the dehydration efficiency of cirrus in the 
tropical tropopause, hence the concentration of water vapor entering the stratosphere*. Cirrus modification has been suggested 
as a geoengineering strategy to counterbalance greenhouse warming’. However, modeling studies show discrepancy on the 
predicted response of cloud properties to the emission of distinct aerosol species” !®-!!, introducing uncertainty in the reliability 
of cirrus geoengineering on climate!”. 

The representation of cirrus in general circulation models (GCMs) remains challenging. Wind velocity, temperature, 
relative humidity and the physicochemical properties of the aerosol may vary at smaller scales (~ 1 to 20 km) than the typical 
resolution used in GCMs (~ 100 — 200 km)!*:!*. Moreover, the ice-nucleating efficiency of aerosol depends on particle size 
and chemistry, and is not completely understood!>. Ice formation is also realized by two different processes: homogeneous ice 
nucleation (HOM), i.e., the spontaneous freezing of supercooled liquid droplets (typically sulfuric acid and sulfate solutions), 
and, heterogeneous ice nucleation (HET), which requires the presence of ice nucleating particles (INP). Although sparse, INP 
can profoundly alter the evolution of clouds®. Much investigation has been devoted over the last decade to elucidating the 
nature of INP showing that a small fraction of atmospheric aerosol (typically some organics, mineral dust and black carbon) are 
capable of nucleating ice'*. The HOM and HET mechanisms can occur simultaneously, and in fact they “compete” during 
cloud formation. HOM nucleation requires higher relative humidity (RH) than HET. Thus, by depleting the increase in RH from 
cooling driven by vertical motion, ice crystals formed by the HET mechanism may inhibit and even prevent HOM nucleation!®. 
However, compared to HOM, HET nucleation leads to lower ice crystal concentration!” !®. 

Vertical motion determines the maximum relative humidity in a cloudy parcel and drives ice nucleation!’. Early studies 
showed that introducing vertical velocity perturbations within parcel model simulations resulted in high variability in ice 
crystal concentration, N;,°?! which has been confirmed by aircraft observations within cirrus*”~**. Modeling studies have 
also shown that low T and high W favors the HOM over the HET mechanism*>-”°, which is reflected in the global distribution 
of N;, simulated in GCMs>!!-*°, Field measurements however suggest predominance of HET and lower N; than implied by 
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Figure 1. Vertical velocity (m s~!) distribution within cirrus at the SGP site from cloud Doppler radar measurements!? (red 
lines) and high resolution simulations (blue bars). Measurements correspond to the period 2000 — 2010 and model results to 
2005 — 2007. Scaling was applied to the model results as described in Methods. Black lines represent the “unscaled 
distribution” from the GSNR. Positive W values indicate updraft. The legends correspond the mean value of 0,, for each case. 


HOM!”:'4-24 consistent with sustained levels of supersaturation found in cirrus clouds”*:*!-3. This suggests that the frequency 
of HOM in cirrus is likely overestimated by GCMs. Recent studies point at overestimation in the parameterized cloud-scale W 
as the likely cause of the discrepancy between models and observations**~*°. 

Because of the separation between the relevant scale for ice nucleation and the scale resolved by the GCMs, it is likely 
that several cloud formation events occur within a model grid cell. This unresolved variability is characterized by a “subgrid” 
distribution of vertical velocity, P(W, Ow), largely determined by its standard deviation, 0,,, since the mean resolved vertical 
transport is slow compared to cloud-scale motions!*:'+79, Cloud-scale motions are typically not resolved by CGMs due 
to their coarse resolution (~ 100 km). To represent subgrid W variability, most GCMs rely on either poorly constrained 
parameterizations or empirical correlations representing particular cloud realizations. A recent study*’” suggested that vertical 
wind velocity may be responsible for about 90% variation in calculated ice crystal formation rates, although it is not clear 
whether the same relation applies to other cloud properties. Field campaign analyses and cloud modeling studies”**-** also 
suggest a strong relation between the effect of aerosol emissions on cloud properties and W. These highlight the importance of 
improving the representation of subgrid W variability in GCMs. 

In this work we develop a method, using ultra high resolution global simulations to directly calculate the global distribution 
of subgrid vertical velocity affecting cirrus formation. By implementing our estimates in a global model and running experiments 
constrained by observations we show that the global distribution of o,, largely determines the balance between homogeneous 
and heterogeneous ice nucleation during the formation of cirrus. 


Methods 


Vertical velocity distribution 
We developed a method to calculate ©(W,o,,) using results from a global climate simulation carried out using the non- 
hydrostatic version of the NASA Goddard Earth Observing System (GEOS-5) over the period 2005-2007°°+!. The simulation 
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Figure 2. Similar to Fig. 1 but for the Manus site. 
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was completed using a cubed-sphere grid with a nominal spatial resolution of 7 km and 72 vertical levels, extending from the 
surface to 0.1 mbar. The time step for integration was set to five minutes and three-dimensional instantaneous meteorological 
fields were saved every hour. The vertical resolution of the model in the upper troposphere is about 0.5 km. This simulation is 
referred as the GEOS-5 “nature run” (GSNR). The output of the simulation amounts to about four petabytes and is used to 
perform Observation System Simulation Experiments*?. The GSNR is capable of resolving mesoscale systems and organized 
convection within large scale midlatitude cyclones, important for the simulation of cirrus*”. 

The global distribution of vertical velocity was calculated collocating W from the 7 km GSNR output to 1° (~ 100 km) 
horizontal resolution (except for the validation studies, for which 0.5° sections were used) so that at least 256 W values from 
the GSNR were used to calculate o,, for each 1° grid cell. This procedure resulted in a hourly, three-dimensional global 
characterization of ®(W,6,,) at the 1°(~ 100km) resolution for the two-year period of the simulation. Monthly averages were 
calculated by averaging 62 using hourly output. Sensitivity tests were carried calculating o,, for 200 km global resolution, 
producing essentially similar results as our 100 km calculation, indicating that very little of the total W variance is resolved at 
scales greater than 100 km. A second sensitivity test was performed collocating results from a short term 3.5 km resolution run 
spanning 15 days in May 2006 to the 1° spatial resolution. After scaling was applied (@ 9 = 1.26, as explained below, Eq. 4), 
this test produced similar o,, as when the 7 km simulation was used (supplemental Figure $1). 

The GSNR resolves waves and vertical motion with periods from minutes to a few hours which are the main drivers of 
cirrus formation?*. However, high frequency waves with periods of a few minutes are responsible for the formation of “pockets” 
of high N; within clouds'+3. The scale of such waves would likely be smaller than the 7 km resolution of the GSNR. To 
account for such a unresolved variability a “scaling” method is developed, as follows. The total spatial vertical velocity variance 
at the 100 km resolution, a. is represented as the sum of the resolved and unresolved components of the high resolution run, 


2 al 2 
On = Ow 7km z. Owunres> (1) 


where en is the spatial variance in W calculated from the 7 km resolution output, and OG, cnies the variance resulting from 
vertical motion with characteristic scales below 7 km. The contribution to the W spatial variability from motion with scales 
greater than 100 km is neglected. This is justified because motion with scales greater than 100 km is resolved at the low 
resolution. To estimate 0} ynres, the approach of Pauluis, et.al. (2006)** is employed. Using a discretized version of the 
equations of motion of an anelastic fluid, the authors derived a relation for the vertical velocity resolved at two different 
horizontal resolutions, e.g., ro and rj, 


Wro I+ RZ = 

where AZ = 6 km corresponds to the resolution at which the increase in kinetic energy from buoyancy is equally distributed 
between its horizontal and vertical components. Pauluis, et.al. (2006)*? showed that Eq. (2) is accurate for horizontal resolutions 
between 2 and 20 km, hence it is suitable to scale the GSNR, 7 km results, to finer resolutions, hence to estimate the variance 
unresolved by the GSNR. In this study we assume rp = 0.1 km as the horizontal scale below which the cloud properties can be 
considered uniform. Using r; = 7 km (..e., the resolution of the GSNR) results in 9 = 1.46. For r; = 3.5 km, Q) 9 = 1.26. 
The choice of rp is rather arbitrary, however for small enough values has little effect on 9 (€.g., @,9 = 1.41 for 79 = 0.5 km). 
Here ro is selected much smaller than the typical cloud size to account for short-lived fluctuations that may affect the relaxation 
of supersaturation in the cloud!+?9:44-49, 

For the special case of a normal distribution at resolution r1, 


W,, i &d (W, Ow.7km) (3) 


where W is the grid-scale vertical velocity. Combining Eq. (2) and Eq. (3) it can be readily shown that, 


Wr ww ® (a1,0W, 1 0 Ow,7km) : (4) 


Which implies 6,, = 01 ,00y,7km. Equation (4) provides a simple argument to scale the resolved variance calculated directly 
from the GSNR output, Oo gins to obtain the total variance driving cloud formation, o, at the 100 km resolution. 

Validation 

The vertical velocity fields from the 7 km simulation were validated against ground-based measurements within cirrus for 
two different sites with diverse topography and synoptic conditions, and for which long term radar retrievals are available!>. 


4/14 


Obs SGP £ 
—@ Model E. 
0.25 
— 0.20 Se 
"e f Si E "e 
& 0.15 E = 
3 E Es 
© 0.10 muege FO 
0.05 
0.00 = 0.00 
Jan'FebMar AprMayJun! Jul IaugSeplOctNovded Jan'FebMarAprMayJun! Jul IsugSeploctNovded 


Figure 3. Monthly mean standard deviation in vertical velocity (m s~!) at the SGP and Manus sites calculated from 7 km 
global output (model) and from radar retrievals (obs). Error bars have been omitted for clarity however the standard error in the 
observations!? is about 0.15 ms~!. 


These correspond to the Southern Great Plains of North America (SGP, 36° 36’ 18” N, 97° 29’ 6” W) and the Manus island 
in the tropical western Pacific (Manus, 2° 3’ 39.64” S, 147° 25’ 31.43” E). SGP is a mid-latitude continental site with large 
variability in temperature and cloud occurrence. Manus is an oceanic site off the coast of Australia with frequent tropical 
convection. Ground-based radar retrievals!* of vertical velocity at each site over the period (2000 — 2010) were obtained from 
the Atmospheric Radiation Measurement program (www.arm.gov/sites). The retrieval algorithm is based on a decomposition of 
the Doppler vertical velocity. The typical error in vertical velocity is about 15 cm s~! for a minimum reflectivity of about —40 
dBz!?. The whole data set spans about 10 years for each site. Data obtained at 10 s intervals for each month were averaged 
over 5 min to match the time step of the GSNR simulation. 

To generate GSNR vertical velocity distributions at the SGP (Fig. 1) and Manus (Fig. 2) sites, instantaneous W values over 
a0.5° x 0.5° area centered at each site were obtained from the 7 km simulation at three hour intervals. This corresponds to 
about 15360 values for each monthly distribution at each site. The model results were restricted to ice mixing ratios above 
5 x 10-> kg kg~! and corresponding ice water content of about 50 mg m~?, selected to match the maximum sensitivity of the 
retrieval method!?. Notice that this includes the spatial and temporal components of the variance since it is difficult to separate 
them in the observations. 

To account for vertical motion with scales below 7 km, likely resulting from in situ turbulence and high-frequency gravity 
waves!*, the distributions were scaled to 0.1 km horizontal resolution using the method outlined in the previous section. 
Equation (2) was applied directly to the simulated distributions, i.e., without invoking the assumption of a normal ® (W, o,,). 
However Figs. 1 and 2 suggest that &(W, o,,) can be adequately approximated using Gaussian functions, in line with literature 
reports2©294648 | 

Comparison of the measured and simulated ®(W,o,,) (Figs. 1 and 2) shows that the G5SNR is capable of generating 
realistic vertical velocity distributions, and in reasonable agreement with observations. At the SGP site it is evident that using 
only the raw GSNR data would result in a too-narrow ® (W, o,,) compared to the measurements (Fig. |, black lines). Scaling 
brings ®(W, o,,) closer to the observations suggesting that a significant fraction of the W variance lies in the sub-7 km range. 
On the other hand, for the Manus site the raw GSNR distribution (Fig. 2, black lines) approximates the observed distribution 
(Fig. 2, red lines), indicating that most of the W variance is resolved at the 7 km resolution. In the latter, scaling may lead to 
overestimation in 6,, since Eq. (2) implicitly assumes that a significant fraction of the W variance lies at the small scales. Thus, 
o,, can be considered a upper limit to the vertical velocity variance. 

Seasonal variation in the large scale environment may lead to differences between the simulated and observed distributions. 
Such deviations are typically within the margin of error of the observations. However they can also signal systematic errors in 
the simulated ® (W, Oy). To study the latter, o,, was calculated for each month over the whole observational record (resulting 
in ten data points per month) and from the GSNR (which are available for two separate years). The results are plotted in Figure 
3. At both sites, the simulated and observed 6,, show little interannual variability, i.e., for the same month the spread in o,, 
between years is typically below 0.05 ms~!. ©, at the SGP site shows a strong annual cycle, whereas at the Manus site it is 
relatively constant over the year. This suggests that location plays an important role in determining o,,. The distinctive behavior 
of 0,, at each site is well represented by the simulation. However the GSNR tends to predict a stronger annual cycle of Oy, at 
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Figure 4. Annual mean standard deviation in vertical velocity (m s~') at a horizontal resolution of 1°, calculated from 7 km 
global output. Left panel: o,, at 250 hPa. Right panel: zonal mean o,,. Maps generated using the NCAR Command Language 
(Version 6.3.0) Software. (2016). Boulder, Colorado: UCAR/NCAR/CISL/TDD. 

http: //dsx.doi.org/10.5065/D6WD3XH5. 


the SGP site than implied by observations, with the maximum o,, occurring too early during the year and an underestimation in 
6,, between August and October, likely related to the low frequency of convective events during the Fall season in the GSNR*°. 


GCM implementation 

The calculated 6,, was used to drive ice nucleation in GEOS-5 and study the impact on the balance between HOM and HET 
processes. Due to the high computational expense of the 7 km simulations, aerosol-cloud interactions are analyzed in the lower 
resolution, 100 km, simulation, but using 0,, obtained from the 7 km run. The goal of the GCM implementation is to analyze 
the statistics of the ice crystal concentration when the GSNR-generated distribution of Oy, is used. The main premises behind 
this approach are that the interannual variability in 0,, is small and that o,, is mostly influenced by local features and convection. 
This is supported by the low interannual variability in o,, found in our 7 km simulation, and in the 10-year-long time series of 
radar retrievals at the SGP and Manus sites (Fig. 3). 

The model’s cloud microphysics and ice nucleation schemes are described elsewhere*’. Briefly, HOM occurs on sulfate 
whereas mineral dust, black carbon and glassy organics are considered INP°’. The ice nucleation rate is weighted by the 
distribution of vertical velocity and explicitly integrated. Large scale deposition, ice cloud fraction, and ice nucleation are 
coupled*’. To minimize the effect of model uncertainties from the transport of aerosol and the meteorological conditions, the 
100 km simulation is constrained using the horizontal wind velocity, temperature and water vapor from the second version of 
the Modern-Era Retrospective Analysis for Research and Applications (MERRA-2)°*”>!, using a relaxation time scale of 6 h. 
Monthly-averaged o,, derived from the GSNR was used to drive ice nucleation. To perform ice nucleation studies, the model 
was run for five years (2005-2010) using a c90 cubed-sphere grid (spatial resolution of about 1°) and 72 vertical levels. Cloud 
physical and optical properties obtained in similar simulations have been shown to be in agreement with available in situ and 
satellite observations*?. 

A cirrus formation event was considered dominated by homogeneous freezing when at least 80% of the ice crystals 
were produced by the HOM mechanism. Modeling studies!’:!® show that cloud formation becomes HET dominated within 
a relative narrow range of INP concentration, typically less than a factor of two. Thus the 80% limit represents the INP 
concentration at which N; becomes strongly affected by HET nucleation. Results using different thresholds, i.e., 20%, 50% 
and 90% are discussed below. Bivariate N; vs. T distributions were calculated by counting the number of data points falling 
within a particular Nj and T combination, using 80 logarithmic bins for N; and 80 linear bins for T. The frequencies were then 
normalized by the most frequent N; vs. T combination within the entire domain. Using this method the expected value of Nj; is 
located around the maximum frequency at each temperature. 
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Figure 5. Probability distribution of o,, for different regions, and for grid cells with no underlying convection. 


Results 


Vertical velocity distribution 

Figure 4 shows the annual global distribution of 0, derived from the GSNR. High values of o,, are found around the Inter- 
Tropical Convergence Zone (ITCZ) and over continental mountain ranges, confirming the notion that underlying convection 
and orography are the main drivers of dynamic variability in the upper troposphere*®-*?->?, For the same reason Oy, is typically 
higher in the Subtropical Northern hemisphere (NH) than in the Southern hemisphere (SH), except over the Andes mountains 
where o,, tends to be high due to orographic uplift. Within the troposphere oy, decreases slightly with height, and becomes 
small above the tropopause due to gravity wave breaking. This is in agreement with observations*> however differs from 
previous work where W, instead of o,,, was assumed to decrease with height"**>. Figure 4 suggests that high values of W do 
occur at high altitude, but they become less frequent near the tropopause since ® (W, 0) becomes narrow. 

Aside from the poles, subtropical regions, usually associated with persistent stratocumulus decks and large scale subsidence 
(i.e., the coasts of Peru, California and Namibia), display the smallest o,, over the globe. These areas are also associated with a 
low frequency of cirrus clouds!. Low oy results in low N; and large ice settling velocities, which may contribute to a more 
efficient cloud removal and explain in part the low cirrus frequency of these regions. At constant pressure 0,, tends to decrease 
poleward from the Tropics in both hemispheres due in part to the breaking of gravity waves at higher pressure levels. In NH the 
lowest 6,, values are found at the northernmost latitudes (above 80°). In SH the lowest values of 6, are found around —60°, 
however 0,, increases south of —60° revealing the orographic effect of the Antarctic continent on W, which may impact the 
formation of polar stratospheric clouds. 

Figure 5 shows the global PDF of o,,, P(o,,), within cirrus, i.e., positive ice water content. As expected, P(o,,) decreases 
exponentially with o,,. However several features of P(o,,) stand out. P(o,) peaks at around 0.02 m s~! and decreases rapidly 
with increasing 6, so that about 90% of the values of o,, are below 0.1 ms~!. This suggests that in most cases cirrus formation 
is driven by large scale vertical transport and inertial gravity waves. Higher values of 6, (up to 0.8 m s~!) are associated with 
high frequency gravity waves from convective and orographic sources*°. Although likely, they are progressively less frequent 
as Ow increases. Globally about one in 10* non-convective cirrus formation events are forced by vertical motion with oy, > 0.8 
ms~!. High values of 6,, are more likely in the Tropics driven by underlying convection. Removing grid-cells with vigorous 
convection results in about a factor of five lower probability of finding o,, > 0.2 ms~!. Values of 0, greater than 0.5 m s~! 
seem to be equally likely in SH and NH indicating that they may result from strong orographic uplift in the mountain ranges. 

There is a slight seasonal variation in the global distribution of o,,. Over the Tropics the highest 0,, coincides with the 
displacement of the ITCZ (Fig. 6). Over the continents and away from orographic features, 6,, is influenced by large scale 
meteorological patterns, i.e., cold fronts, midlatitude cyclones, and subtropical jets?’. In the subtropics 0,, tends to be larger 
during the summer months, particularly over land. The two years of the GSNR simulation show similar patterns: highest in the 
ITCZ and in orographic regions and lowest in the high latitudes with little variation in the Subtropical regions. Location seems 
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Figure 6. Difference between the monthly and annual o,, at 100 km resolution calculated from the GSNR 7 km resolution 


output at the 250 hPa pressure level. Maps generated using the NCAR Command Language (Version 6.3.0) Software. (2016). 
Boulder, Colorado: UCAR/NCAR/CISL/TDD. http: //dx.doi.org/10.5065/D6WD3XH5. 
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Figure 7. Global distribution of the frequency of cirrus events dominated by homogeneous ice nucleation, vertically weighted 
by cloud fraction (a) and zonal mean (b). Maps generated using the NCAR Command Language (Version 6.3.0) Software. 
(2016). Boulder, Colorado: UCAR/NCAR/CISL/TDD. http: //dx.doi.org/10.5065/D6WD3XH5. 


to be the main factor determining 0,,. This is supported by the data at the Manus and SGP sites (Fig. 3), and it is in line with 
literature reports showing low interannual variability in ®(W, o,,) for the same location**:>’. A longer simulation period is 
however required to further study the interannual variability in ® (W,0,,) and will be the subject of a future study. 

The results shown in Fig. 4 can be used to evaluate the parameterization of o,, currently used in GEOS-5 (Barahona et al. 
2014”, c.f. Fig. 4). Compared to the parameterized results, 6,, is lower in the marine midlatitudes and higher in the Tropics. 
This parameterization used the orographic stress and turbulence to estimate 0,,. Over the ocean the total surface stress was 
used, leading to an overestimation of Oy in the midlatitudes. On the other hand, the parameterization did not explicitly account 
for gravity waves generated by convection leading to underestimation of o,, in the Tropical region. Over land the agreement is 
better since orographic features drive gravity wave generation>”. As similar parameterizations are used in many GCMs, this 
exercise shows that Fig. 4 may provide a way to evaluate their performance. 


Implications for ice cloud formation 

The annual mean frequency of homogeneous freezing events during cirrus formation calculated with GEOS-5 (forcing ice 
nucleation with oy calculated from the GSNR output) is shown in Figure 7. High o,, and low T tend to lead to high annual 
HOM frequency. This is because under such conditions a larger INP concentration is required to offset the increase in RH from 
expansion cooling, and is evidenced by the correspondence in the spatial patterns seen in Figs. 4 and 7. In the Tropical region 
where high values of 6,, are more likely, HOM dominates about 50 — 60% of the cloud formation events, and up to 80% in 
the coldest regions of the Tropical tropopause, over the Central Pacific, where a lack of INP also contributes to diminish the 
frequency of HET nucleation. Outside the Tropical and Subtropical regions, the frequency of cloud events dominated by HOM 
nucleation is generally low (below 30%), particularly in the NH. In fact, in the Arctic clouds tend to form almost exclusively 
by HET (HOM frequency is below 10%), mostly resulting from the low o,, which precludes efficient HOM nucleation, even 
though the INP availability is low. HOM-dominated cirrus events are also less frequent over North Africa due to the presence 
of dust and black carbon, and to the absence of orographic features and convection that would produce high o,,. 

While the global mean HOM-dominated frequency is relatively constant over the year at ~ 37%, in both the SH and NH 
there are peaks in the winter months, since low temperature favors HOM (Fig. 8). Maximum HOM occurrence in the SH is 
45% during the winter months, while it is only 30% in the NH winter. Even though o,, is higher on average in the NH than in 
the SH, HOM is more prevalent in the latter during most of the year since INP like mineral dust are abundant in the NH°®. 
The seasonal differences in the SH are more pronounced than in the NH due to the larger temperature and aerosol variability, 
and lower 6, in the former. In the Tropics, a seasonal cycle is also present reflecting the strength and position of the ITCZ 
influencing o,, and the annual variation in black carbon and dust concentration. 

The balance between HOM and HET during cloud formation is significantly influenced by the concentration of ice 
nucleating particles, Nijnp. Supplementary Figure S2 shows the frequency of Mjnp calculated at the maximum RH during 
cloud formation, RHmax. Ninp increases steeply around RHmax = 15% to a maximum global average value of about 3 L-! at 
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Figure 8. Monthly mean homogeneous ice nucleation frequency for the Tropical (latitude —30° to 30°) and the Northern (NH, 
latitude 30° to 60°), and Southern (SH, latitude —30° to —60°) extratropical regions. 


RH max = 30% (~ 5 L~! for NH). The apparent decrease in Mnp at high RHmax is caused by competition between HOM and 
HET nucleation. For RA;yq, < 40% where HOM does not occur, the simulated N;yp shows similar characteristics as available 
reports>*. This is in part by design as the heterogeneous ice nucleation spectrum used in the model is derived from a collection 
of available field data>®. Figure $2 shows that GEOS-5 simulates realistic Ninp statistics. 

Globally, about 70% of the time cirrus form in situations where HOM and HET occur simultaneously, with HET being 
more prevalent. This finding is at odds with previous modeling studies where HOM is the predominant cirrus formation 
mechanism*-*”:>*. Our results are however in better agreement with field campaign data suggesting a significant role of dust 
and other INP species in cirrus formation!”, and, lower N; and higher saturation ratios than implied by HOM®:'*. This suggests 
that the parameterization of o,, may be responsible for the higher HOM frequency typically found in modeling studies. Notably, 
forcing ice nucleation with our estimate of 0, also results in good agreement of N; with field campaign data at very low 
temperature (Fig. 9, T < 200 K), where most modeling studies show high overestimation of N,>*°. This suggests that the 
overestimation in many models may be a result of poorly constrained 0,,. GEOS-5 however seems to overestimate the frequency 
of low Nj; for T > 230 K. It must be noticed that ice shattering on in-situ probes, leading to overestimation in the in-situ N;, are 
a likely artifact of the measurements at such temperatures®. Results for the NH, where most cirrus field campaigns have taken 
place** show similar tendencies (see supplementary Figure S3) with slightly lower variability in N; and better agreement with 
observations for T > 230K 

The effect of selecting different thresholds to define a HOM-dominated cirrus is shown in the supplementary Figure S4. 
Changing the definition of a HOM-dominated cloud from 90% to 50% of N; produced by HOM, increases the HOM frequency 
from 33% to 43%. The latter corresponds to the minimum threshold where a cloud could be referred as HOM-dominated since 
for any value below 50% most of the ice crystals would in fact be produced by HET nucleation. Even at the 50% threshold, 
the HOM frequency is still much lower than reported values showing that our conclusions are not significantly influenced by 
the selected threshold. Using a more strict definition, where HOM still occurs but is not dominant, i.e., 20% of ice crystals 
produced by HOM (actually a HET-dominated cloud) leads to global HOM frequency of about 48% which is still low compared 
to current reports” *?, 


Discussion and Conclusions 


This work reports for the first time the direct estimation of the subgrid spatial variance in vertical velocity at the scale relevant 
for cirrus formation. Figure 4 shows that 0, has considerable spatial and seasonal variability. It is important for GCM 
parameterizations of o,, to reproduce such a spatial variability. This indicates that parameterizations based on individual field 
campaigns should be used with caution when applied to the global scale. Large-scale dynamics, turbulence, orographic uplift, 
and underlying convection impact o,,. Even State-of-the-Art parameterizations of o,*°'>* neglect the effect of convection 
generating gravity waves that can impact cirrus formation, which may result in a too-low frequency of high o,, as indicated by 
Fig. 5. 

Using the direct estimate of oy, to drive the GEOS-S ice nucleation scheme results in a lower predominance of homogeneous 
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Figure 9. Global frequency distribution of in-cloud ice crystal number concentration as a function of temperature from 
GEOS-5 output over a 2-year subset (2005-2006). Solid lines represent the 25 % and 75 % quantiles from a compilation of 


field campaign observations». 


ice nucleation than previously simulated, in better agreement with field measurements. Such lower frequency of HOM events 
also results in a better agreement of simulated N; with observations, particularly at low T (Fig. 9). This is significant as cloud 
formation theory typically predicts high N; for low temperature*”’’, a trend also found in GCM studies. One way to reconcile 
theory and observations is to assume very low W at low T°°-°®>4, which however conflicts with observations of high W (> 0.5 
m s~!) near the tropopause!*?-°!, Our results provide a way to reconcile these two views. High values of W do occur at low 
T, however the low 0, at high altitude limits the frequency of cloud formation events driven by high W, and on average leads 
to lower N;. This is agreement with previous work suggesting the structuring of cirrus by dynamics” **-*> and an episodic 
nature of HOM in cirrus!*. 

Our new estimate provides a way to validate new parameterizations of 0,, at a global scale. Several uncertainties however 
remain in the modeling of cirrus clouds, the most significant being the concentration of INP in the upper troposphere. Even 
though our results are in relative good agreement with available reports (Supplementary Figure S2), still few studies provide 
observations of Njyp with a wide range of aerosol concentrations, temperature and relative humidity and enough spatial 
coverage that would allow comprehensive validation of GCM predictions. The characterization of W at low T also requires 
better techniques with smaller errors, since in some cases they can be as high as o,,'°. In turn an improved representation of 
0,, in GCMs may help reducing the uncertainty surrounding the estimation of the impact of aerosol emissions on cirrus, and 
eventually lead to a better prediction of future climate. 
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Figure S1. Monthly mean standard deviation in vertical velocity (m s~!) at pressure level 250 hPa for May 2006 at horizontal 
resolutions of 100 km and 200 km. The top panel was calculated using the 7 km global output. The middle and bottom panels 
were calculated using 3.5 km global output. Maps generated using the NCAR Command Language (Version 6.3.0) Software. 


(2016). Boulder, Colorado: UCAR/NCAR/CISL/TDD. http: //dx.doi.org/10.5065/D6WD3XH5. 
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Figure S2. Frequency distribution of the concentration ice nucleating particles, Njnp, as a function of the maximum in-cloud 
relative humidity with respect to ice. GEOS-5 output was used over a 2-year subset (2005-2006) and 100 km horizontal 
resolution. Left panel: Global. Right panel: Northern Hemisphere. 
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Figure S3. Frequency distribution of in-cloud ice crystal number concentration as a function of temperature for the Northern 
Hemisphere, from GEOS-5 output over a 2-year subset (2005-2006). Solid lines represent the 25 % and 75 % quantiles from a 
compilation of field campaign observations!. 
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Figure S4. Global distribution of the frequency of cirrus events with HOM nucleation producing 90% (a, b), 50% (c, d), and 
20% (e, f) of ice crystals. Left panels: vertically weighted by cloud fraction. Right panels: zonal mean. Maps generated using 
the NCAR Command Language (Version 6.3.0) Software. (2016). Boulder, Colorado: UCAR/NCAR/CISL/TDD. 
http://dx.doi.org/10.5065/D6WD3XH5. 
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